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Abstract 

The extended Airy kernel describes the space-time correlation functions for the 
Airy process, which is the hmiting process for a polynuclear growth model. The 
Airy functions themselves are given by integrals in which the exponents have a 
cubic singularity, arising from the coalescence of two saddle points in an asymp- 
totic analysis. Pearcey functions are given by integrals in which the exponents 
have a quartic singularity, arising from the coalescence of three saddle points. A 
corresponding Pearcey kernel appears in a random matrix model and a Brownian 
motion model for a fixed time. This paper derives an extended Pearcey kernel by 
scaling the Brownian motion model at several times, and a system of partial differ- 
ential equations whose solution determines associated distribution functions. We 
expect there to be a limiting nonstationary process consisting of infinitely many 
paths, which we call the Pearcey process, whose space-time correlation functions 
are expressible in terms of this extended kernel. 

I. Introduction 

Determinantal processes are at the center of some recent remarkable developments in 
probability theory. These processes describe the mathematical structure underpinning 
random matrix theory, shape fluctuations of random Young tableaux, and certain 1 + 1 
dimensional random growth models. (See [21 [HI 1101 HH 120] for recent reviews.) Each 
such process has an associated kernel K{x,y), and certain distribution functions for 
the process are expressed in terms of determinants involving this kernel. (They can 
be ordinary determinants or operator determinants associated with the corresponding 
operator K on an space.) Typically these models have a parameter n which might 
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measure the size of the system and one is usually interested in the existence of limiting 
distributions as —>■ oo. Limit laws then come down to proving that the operator 
Kn, where we now make the n dependence explicit, converges in trace class norm to 
a limiting operator K. In this context universality theorems become statements that 
certain canonical operators K are the limits for a wide variety of Kn. What canonical 
K can we expect to encounter? 

In various examples the kernel Kn{x, y) (or, in the case of matrix kernels, the matrix 
entries Kn,ij{x,y)) can be expressed as an integral 

To study the asymptotics of such integrals one turns to a saddle point analysis. Typically 
one finds a nontrivial limit law when there is a coalescence of saddle points. The simplest 
example is the coalescence of two saddle points. This leads to the fold singularity 
02 (-2) = + \z m. the theory of Thom and Arnold and a limiting kernel, the Airy 
kernel or the more general matrix- valued extended Airy kernel ^3 • 

After the fold singularity comes the cusp singularity 03(2;) = + X^z^ + \\z. 
The diffraction integrals, which are Airy functions in the case of a fold singularity, now 
become Pearcey functions ^6j. What may be called the Pearcey kernel, since it is 
expressed in terms of Pearcey functions, arose in the work of Brezin and Hikami [HI U\ 
on the level spacing distribution for random Hermitian matrices in an external field. 
More precisely, let if be an n x n GUE matrix (with n even), suitably scaled, and Hq 
a fixed Hermitian matrix with eigenvalues ±a each of multiplicity n/2. Let n —>■ 00. If 
a is small the density of eigenvalues is supported in the limit on a single interval. If a 
is large then it is supported on two intervals. At the "closing of the gap" the limiting 
eigenvalue distribution is described by the Pearcey kernel. 

Bleher, Kuijlaars and Aptekarev [H 0] have shown that the same kernel arises in 
a Brownian motion model. Okounkov and Reshetikhin J3] have encountered the same 
kernel in a certain growth model. 

Our starting point is with the work of Aptekarev, Bleher and Kuijlaars With 
n even again, consider n nonintersecting Brownian paths starting at position at time 
r = 0, with half the paths conditioned to end at 6 > at time r = 1 and the other 
half conditioned to end at —b. At any fixed time this model is equivalent to the random 
matrix model of Brezin and Hikami since they are described by the same distribution 
function. If b is of the order n^/^ there is a critical time such that the limiting 
distribution of the Brownian paths as n — 00 is supported by one interval for r < 
and by two intervals when r > t^. The limiting distribution at the critical time is 
described by the Pearcey kernel. 

It is in searching for the limiting joint distribution at several times that an extended 
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Pearcey kernel arises.^ Consider times < ri < ■ • • < < 1 and ask for the probability 
that for each k no path passes through a set at time r^. We show that this probabihty 
is given by the operator determinant det(J — Kx) with an m x m matrix kernel K{x, y), 
where X{y) = diag (Xx^d/))- 

We then take b = n^l"^ and scale all the times near the critical time by the substitu- 
tions Tfc 1/2 + n'^^'^Tk and scale the kernel by x — > n~^^^x, y n~^^'^y. (Actually 
there are some awkward coefficients involving 2^/^ which we need not write down ex- 
actly.) The resulting limiting kernel, the extended Pearcey kernel, has i,j entry 

3_ f g-sV4+T,s2/2-ys+tV4-T,i2/2+a;t dt 

47r^ Jc J-ioo S — t 

plus a Gaussian when i < j. The t contour C consists of the rays from ±ooe*'^/^ to 
and the rays from to ±e~*'^/^. For m = 1 and ri = this reduces to the Pearcey kernel 
of Brezin and Hikami.^ 

These authors also asked the question whether modifications of their matrix model 
could lead to kernels involving higher-order singularities. They found that this was so, 
but that the eigenvalues of the deterministic matrix Hq had to be complex. Of course 
there are no such matrices, but the kernels describing the distribution of eigenvalues 
oi Hq + H make perfectly good sense. So in a way this was a fictitious random matrix 
model. In Section V we shall show how to derive analogous extended kernels and limiting 
processes from fictitious Brownian motion models, in which the end-points of the paths 
are complex numbers. 

For the extended Airy kernel the authors in (21] derived a system of partial differential 
equations, with the end-points of the intervals of the Xk as independent variables, whose 
solution determines det(J — Kx).^ Here it is assumed that each Xk is a finite union of 
intervals. For m = 1 and Xi = (^, oo) these partial differential equations reduce to 
ordinary differential equations which in turn can be reduced to the familiar Painleve II 
equation. In Section IV of this paper we find the analogous system of partial differential 
equations where now the underlying kernel is the extended Pearcey kernel.'^ Unlike the 
case of the extended Airy kernel, here it is not until a computation at the very end that 
one sees that the equations close. It is fair to say that we do not really understand, from 
this point of view, why there should be such a system of equations. 

^It was in this context that the extended Airy kernel, and other extended kernels considered in |21j . 
arose. 

^In the external source random matrix model, an interpretation is also given for the coefficients of 

and in the exponential. It is not related to time as it is here. 
Equations of a different kind in the case m = 2 were found by Adler and van Moerbeke [T]. 

^In the case m = 1 the kernel is integrable, i.e., it is a finite-rank kernel divided hy x — y. (See 
footnote [7| for the exact formula.) A system of associated PDEs in this case was found in [7j, in the 
spirit of |19). when Xi is an interval. This method does not work when to > 1, and our equations are 
completely different. 
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The observant reader will have noticed that so far there has been no mention of the 
Pearcey process, supposedly the subject of the paper. The reason is that the existence of 
an actual limiting process consisting of infinitely many paths, with correlation functions 
and spacing distributions described by the extended Pearcey kernel, is a subtle proba- 
bilitic question which we do not now address. That for each fixed time there is a limiting 
random point field follows from a theorem of Lenard [1^1 IT^ (see also PH]); since that 
requires only a family of inequalities for the correlation functions which are preserved 
in the limit. But the construction of a process, a time-dependent random point field, is 
another matter. Of course we expect there to be one. 

II. Extended kernel for the Brownian motion model 

Suppose we have n nonintersecting Brownian paths. It follows from a theorem of 
Karlin and McGregor ^21 that the probability density that at times tq, . . . , r^+i their 
positions are in infinitesimal neighborhoods of Xoi, . . . , Xm+i,i is equal to 

m 

\Y det {P{Xm,i, Xm+l,j, am)) (2.1) 
fc=0 

where 
and 

P{x,y,a) = (7ra)-^/2g-(x-j/)V-. 
The indices i and j run from to n — 1, and we take 

To = 0, Tm+1 = 1. 

We set all the xoi = a, and Xm+ij = bj, thus requiring our paths to start at Oi and 
end at bj. (Later we will let all Oj — 0.) 

By the method of jH] (modified and somewhat simplified in |21j) we shall derive an 
"extended kernel" K{x, y), which is a matrix kernel {Kk£{x, y))'^^^^^ such that for general 
functions fi,...,fm the expected value of 

m n—l 

n 11(1 +^(^^^)) 

k=l i=0 

is equal to 

det (I-Kf), 

where 

fiy) = diag (My)). 
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In particular the probability that for each k no path passes through the set at time 
Tfc is equal to 



(The same kernel gives the correlation functions |S] . In particular the probability density 
dTU is equal to {nl)-"" det{Kke{xki, Xij))k,e=i,...,m; i,j=o,...,n-i-) 

The extended kernel K will be a difference H — E, where E is the strictly upper- 
triangular matrix with k,i entry P{x,y, — t^) when k < i, and where Hke{x,y) is 
given at first by the rather abstract formula ()2.5p below and then by the more concrete 
formula ()2.6|) . Then we let all — and find the integral representation ()2.11|) for the 
case when all the Brownian paths start at zero. This representation will enable us to 
take the scaling limit in the next section. 

We now present the derivation of K. Although in the cited references the deter- 
minants at either end (corresponding to A; = and m in 1)2.111 ) were Vandermonde 
determinants, it is straightforward to apply the method to the present case. Therefore, 
rather than go through the full derivation again we just describe how one finds the 
extended kernel. 

For i, j = 0, . . . ,n — l we find Pi{x), which are linear combination of the P{x, ak, (Tq) 
and Qj{y), which are linear combination of the P{y,bk, (Jm) such that 



det (I-KX), 



where 



X(i/) = diag (XxM)- 




m— 1 



k=l 



Because of the semi-group property of the P{x,y, r) this is the same as 




(2.2) 



We next define for k < £ 




r=k 



e-1 = 



P{xk,xe, Te - Tk). 



Set 



Pi — Pli 



Q 



and for A; > 1 define 




(2.3) 
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and for k < m define 

Qkj{x) = I Ekm{x, v) Qj{v) dv = I P{x, V, r„ - r^) Qj{v) dv. (2.4) 



(These fiold also for Pu and Qmj if we set Ekk{x,y) = 6{x — y).) 
The extended kernel is given hj K = H ~ E where 



Hke{x, y) = Y1 Qki{x) Puiy), (2.5) 

1=0 

and Eke{x, y) is as given above for /c < £ and equal to zero otherwise. 

This is essentially the derivation in |8j applied to this special case. We now determine 
Hkii{x, y) explicitly. 

Suppose 

Pi{^) = ^PikP{x,ak, ao), 

k 



If we substitute these into ()2.2|) and use the fact that ctq + — ri + (Jm = 1 we see that 
it becomes 



27r fc^^ 



Thus, if we define matrices P, Q and A by 



t 



then we require PAQ = 
Next we compute 

PrAy) = / P{y, U, Tr - Ti) Pi{u) du = ^Pik P{y, ak, Tr), 
'' k 

Qsj{x) = / P{X, V, Tm - Ts) Qj{v) dv = Qji P{x, h, 1 - T-,). 

Hence 

Hrs{x,y) = J2Qriix) Psi{y) = ^P(a;,6^, 1 - Tr) qapik Piy,ak, r,). 

i i,k,£ 

The internal sum over i is equal to the i, k entry of Q^P = \p2/n A^^ . So the above can 
be written (changing indices) 



If we set B — (e^"'''^ ) then this becomes 

Hke{x, y) = V2^Yl ^(^' 1 - r^) e'" {B-% e^' P{y, a„ r,). (2.6) 



This gives the extended kernel when the Brownian paths start at the aj. Now we are 
going to let all a-, — > 0. 

There is a matrix function D — D{ao, . . . , a„_i) such that for any smooth fuction /, 



lim L'(ao, . . • , a„_i) 



/ /(«o) \ 



V /(«n-l) / 



m 



Here lima.^o is short for a certain sequence of limiting operations. Now B ^ applied to 
the column vector 

{e^"p{y,aj,n)) 
equals {DB)~^ apphed to the vector 

L'(e«?P(y,a,-, T,)). 

When we apply lima._^o this vector becomes 

{die'^'/'P{y,a,Te)\a=o), 

while DB becomes the matrix ((2 6^)'), which is invertible when all the bj are distinct. 
If wc set V = (b/) then the hmiting {DB)~^ is equal to diag (2"-'). Thus we have 
shown that when all Ui — 

Hu{x.y) - a/2^ ^ P{x, bi, 1 - Tk) e'" {V-% di e'^' P{y, a, r,) U=o. (2.7) 
The next step is to write down an integral representation for the last factor. We have 

I _jP_ fioo _[^s^+2s(a--JL-] 



e" P{y, a, n) 



\/2(i - n) 



ds. 



Hence 



m^2{l - re) 



gi-x, / gJei-^^ ds. (2.8) 
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Next we are to multiply this by {V and sum over j. The index j appears in ()2.8p 
only in the factor in the integrand, so what we want to compute is 

Y.iV'%s^. (2.9) 

j=0 

Cramer's rule in this case tells us that the above is equal to Aj/A, where A denotes the 
Vandermonde determinant of 6 = {60, • • • , &n-i} and Aj the Vandermonde determinant 
of b with bi replaced by s. This is equal to 

-pj- s - br 

J-J. bi - br ■ 



Observe that this is the same as the residue 




(2.10) 



This allows us to replace the sum over i in ()2.7|) by an integral over t. In fact, using (j2l 
and the identification of ()2.9|1 with ()2.1()|1 we see that the right side of ()2.7j) is equal to 

1 1 j!_ r , .2 -L^ s^-^^ -n s — br ds dt 



2n Jtt{1- Ti) ■^C J-ioo \^t-br s-t 



where the contour of integration C surrounds all the bi and lies to one side (it doesn't 
matter which) of the s contour. Thus 





1 — Te 1 — Tf. 



V(l - n) (1 - r.) 

-'OC. ^fc ^2 I 2xt I T,, ^2 2sy _^s — br dsdt 



X / e"~'+~+~' -~ TT^-^ (2.11) 

Jc J-ioo r t — br S — t 

In this representation the s contour (which passes to one side of the closed t contour) 
may be replaced by the imaginary axis and C by the contour consisting of the rays from 
iooe*'^/^ to and the rays from to ±ooe~*'^/^. (We temporarily call this the "new" 
contour C.) To see this let Cr denote the new contour C, but with R replacing 00 and 
the ends joined by two vertical lines (where has positive real part). The t contour may 
be replaced by Cr if the s contour passes to the left of it. To show that the s contour 
may be replaced by the imaginary axis it is enough to show that we get when the 
s contour is the interval [—iR, iR] plus a curve from iR to —iR passing around to the 
left of Cr. If we integrate first with respect to s we get a pole at s = t, and the resulting 
t integral is zero because the integrand is analytic inside Cr. So we can replace the s 
contour by the imaginary axis. We then let — >^ cxd to see that Cr may be replaced by 
the new contour C. 
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III. The extended Pearcey kernel 



The case of interest here is where half the br equal b and half equal —b. It is convenient 
to replace n by 2n, so that the product in the integrand in ()2.1H1 is equal to 

/ 2 7 2 \ " 



We take the case b = n^l'^ . We know from P] that the critical time (the time when 
the support of the limiting density changes from one interval to two) is 1/2, and the 
place (where the intervals separate) is x = 0. We make the replacements 

and the scaling 

-1/4 -1/4 

X ^ n 'X, y n ' y. 

More exactly, we define 

Kr.,,{x,y)=n-^'^K,^{n'^/^x, n'^^y), 

with the new definition of the r^. (Notice the change of indices from k and I io i and j. 
This is for later convenience.) 

The kernel En,ij{x,y) is exactly the same as Eij{x,y). As for Hn^ij{x,y), its integral 
representation is obtained from ()2.11|) by the scaling replacements and then by the 
substitutions s — > n^^^ s, t — > n^^^ t in the integral itself. The result is (we apologize for 
its length) 



Hn,ijix,y) = 2 1, 



2y^ 2x^ 



7r2 ^(1 _ 2n-y^Ti) (1 - 2n-y^Tj) I ^'^^ " ^^j - 2n 



lOO 



X / / exp<^-r2"/"- T—r77^r + 



C J —ioo 



1 - 2n-i/2r,- 1 - 2r2-i/V,- 



r 1 + 2n-VVj 2 4t/g ] dsdt 

xexpjn ^ _ 2^-1/2^^. ^ 1-2^-1/2^.^/ \l-t^/n^/^j s-f 

We shall show that this has the limiting form 

1 r „4/o,^^.„2 ,.,„,^4/o A ^.,2 , A dsdt 



(3.1) 



^Pearcey ^ __}_ /" ^-g4/2+4r,- s^-4iL;s+t4/2-4 t^+4a;t ""^ 2) 

' TT^ JC J -ioo S —t^ 
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where, as in fl2.11|) and ()3.1|) . the s integration is along the imaginary axis and the 
contour C consists of the rays from ±ooe^'^^'^ to and the rays from to ±ooe~^'^^^. 
Precisely, we shall show that 

Jirn Hn,ijix, y) = H^P''''^{x, y) (3.3) 

uniformly for x and y in an arbitrary bounded set, and similarly for all partial deriva- 
tives.^ 

The factor outside the integral in 1)3.11) converges to — l/vr^. The first step in prov- 
ing the convergence of the integral in ()3.1|) to that in ()3.2|) is to establish pointwise 
convergence of the integrand. 

The first exponential factor in the integrand in p.ip is 

exp { - ^ ^ 0(^-1/2)) ^2 ^ (4 + 0(^-1/2)) 

while the second exponential factor is 

exp {{n^/^ + 4tj + 0{n-'^^)) - (4 + 0{n-^/^)) ys}. (3.4) 

When s = o{n^^^), t = o(n^/^) the last factor in the integrand is equal to 

exp ^2 ^ ^4/2 + o(tVn) - n'/^ - + o{s^/n)]. 
Thus the entire integrand (aside from the factor l/(s — t)) is 

exp { - (1 + 0(1)) + (4r, + 0{n-^'^)) - (4 + 0{n-^'^)) ys] 

X exp {(1 + 0(1)) tV4 - (4ri + 0{n-^l^)) + (4 + 0{n-^'^)) xt] (3.5) 

In particular this establishes pointwise convergence of the integrands in ()3.H) to that 
in ()3.2|) . For the claimed uniform convergence of the integrals and their derivatives it is 
enough to show that they all converge pointwise and boundedly. To do this we change 
the t contour C by rotating its rays slightly toward the real line. (How much we rotate 
we say below. We can revert to the original contour after taking the limit.) This is 
so that on the modified contour, which we denote by C", we have 3fJt^ > as well as 

m'^ < 0. 

The function l/{s — t) belongs to L'^ for any g < 2 in the neighborhood of s = t = 
on the contours of integration and to L'^ for any q > 2 outside this neighborhood. To 
establish pointwise bounded convergence of the integral it therefore suffices to show that 

^The constants in (|3.2|l are different from tfiose in 1)1. a matter of no importance. In the next 
section we sfiall make the appropriate change so that they agree. 
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for any p G (1, oo) the rest of the integrand (which we know converges pointwise) has 
norm which is uniformly bounded in x and y.^ The rest of the integrand is the product 
of a function of s and a function of t and we show that both of these functions have 
uniformly bounded norms. 

From ()3.5p it follows that for some small e > the function of t is 0(e^*'/2+o(|t|2)) 
uniformly in x if \t\ < ev}l^. Given this e we choose C to consist of the rays of C 
rotated slightly toward the real axis so that if 6* = argt^ when t & C then cos 6' = 6^/2. 

When t e C and |t| < en^'^ the function of t is 0(e^*'/2+o(|t|2)^ _ Q(^^cos2e\t\'^/2y 
Since cos 2^ < the U' norm on this part of C is 0(1). 

When t G C and |t| > en^^^ we have 

|1 - tVnl/2|2 = 1 + ^-l|t|4 _ 2^-1/2 cOs9\t\\ 

But 

n-l|t|^ - 272-1/2 ^^^0 |^|2 ^ ^-1/2|^|2 (^-1/2|^|2 _ ^2) > q 

when \t\ > en^l^. Thus \\-t^ln^l'^\ > 1 and the function of t is 0(6"™'^^'^'^' = 

We have shown that on C the function of t has uniformly bounded norm. 

For the norm of the function of s we see from ()3.4j) that the integral of its pth 
power is at most a constant independent of y times 



JO 



(We replaced s by is, used evenness, and took any r > —ApTj.) The variable change 
s — > n^/^ s replaces this by 

1/4 /""^ g-pn(s2-log(l+s2))+jil/2^s2 



n 

The integral over (1, cxd) is exponentially small. Since — log(l + s^) > s^/2 when 
s < 1, what remains is at most 

^1/4 f g-pnsV2+0(nV2,2) ^ ^ g-p.V2+0(s^) 

which is 0(1). 



®That this suffices follows from the fact, an exercise, that if {fn} is a bounded sequence in L'p 
converging pointwise to / then (/„, g) — > (/, g) for all g S L'^, where p = q/{q — 1). We take /„ to be 
the integrand in except for the factor l/(s — i) and g to be l/(s — i), and apply this twice, with 
q < 2 in a neighborhood of s = t = and with q > 2 outside the neighborhood. 



11 



This completes the demonstration of the bounded pointwise convergence of Hn,ij{x, y) 
to H^^^^^'^^ . Taking any partial derivative just inserts in the integrand a polynomial in 
X, s and t, and the argument for the modified integral is virtually the same. This 
completes the proof of 



IV. Differential equations for the Pearcey process 



We expect the extended Pearcey kernel to characterize the Pearcey process, a point 
process which can be thought of as a infinitely many nonintersecting paths. Given sets 
Xfc, the probability that for each k no path passes through the set Xk at time is equal 
to 

det {I-KX), 

where 

= diag iXxM)- 

The following discussion follows closely that in pP . We take the case where each 
is a finite union of intervals with end-points ^kw) w = 1, 2, . . ., in increasing order. If 
we set 

R={I -KXY^K, 

with kernel R{x,y), then 



dk^det {I-KX) = (-ir+'i?fefc(a^,a«,)- 

(We use the notation dkw for S^j,^.) We shall find a system of PDEs in the variables ^kw 
with the right sides above among the unknown functions. 

In order to have the simplest coefficients later we make the further variable changes 

s ^ s/2'/\ t ^ 

and substitutions 

x^2i/V4, y^2VV4, r, ^ 2^/ VS- 
The resulting rescaled kernels are (we omit the superscripts "Pearcey") 



47r2 

which is and 



Eij{x,y) 



2tT {tj - Ti] 
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Define the vector functions 

We tliink of <y9 as a column m-vector and ■i/' as a row m-vector. Tlieir components are 
Pearcey functions/ The vector functions satisfy the differential equations 

^"\x) -r^\x) + xv9(x) = 0, ^"\y) - ^\y) r -y^{y) = 0, 

where r = diag {jk). 

Define the column vector function Q and row vector function P by 

Q = {I-KX)-^^, P = ij{I -XK)-\ (4.1) 

The unknowns in our equations will be six vector functions indexed by the end-points 
kw of the Xfc and three matrix functions with the same indices. The vector functions 
are denoted by 

q\ q\ p, v\ v" ■ 

The first three are column vectors and the second three are row vectors. They are defined 
by 

and analogously for p, p', p" . The matrix function unknowns are 

1 'XI ' y 

defined by 

(Here Rxij, for example, means dxRij{x,y).) 

The equations themselves will contain the matrix functions rxx, fxy, fyy defined 
analogously, but we shall see that the combinations of them that appear can be expressed 
in terms of the unknown functions. 

The equations will be stated in differential form. We use the notation 
^ = diag i^kw), = diag (d^kw), s = diag ((-1)^"+^). 



^ In case m — 1 the kernel has the exphcit representation 



K{x,y) = 



x-y 



in terms of the Pearcey functions. (Here we set r = ri. This is the same as the i, i entry of the matrix 
kernel if we set r = r^.) This was shown in Another derivation will be given in footnote 1^1 
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Recall that g is a column vector and p a row vector. 
Our equations are 



dr 


= —rsd^r + d^r^ + ry d^, 


(4.2) 


dr^ 


= -r^sd^r + d^r^^ + r^yd^, 


(4.3) 


dry 


= -rsd^ry + d^ r^y + ryy d^, 


(4.4) 


dq 


= dC^ q' — r s d^ q, 


(4.5) 


dp 


= p' dC, — pd^ sr, 


(4.6) 


dq 


= d^ q" -r^sdi q, 


(4.7) 


dp 


= p'd^ — pd^sTy, 


(4.8) 


dq" 


= d^ {t q' - ^ q + r s q" - Ty s q' + Tyy s q - r T s q) - r^^ s d^ q, 


(4.9) 


dp" 


= {p T + p ^ + p" s r — p' s r^ + p s r^x —psTr)d^ — pd^sryy. 


(4.10) 



One remark about the matrix r in equations ()4.9|) and ()4.10|) . Earlier r was the 
mxm diagonal matrix with k diagonal entry r^. In the equations here it is the diagonal 
matrix with kw diagonal entry r^. The exact meaning of r when it appears will be clear 
from the context. 

As in j25) what makes the equations possible is that the operator K has some nice 
commutators. In this case we also need a miracle at the end.^ 

Denote hj ip^ijj the operator with matrix kernel {^Pi{x) tpjly)), where ipi and ipj are 
the components of ip and respectively. 

If we apply the operator dx + dy to the integral defining Hij{x,y) we obtain the 
commutator relation [D, H] = —ip^ip, where D = d/dx. Since also [D, E] = we have 

[D, K] = -ip^tjj. (4.11) 

This is the first commutator. From it follows 

[D, Kx] = -^®i)X + K5. 

Here we have used the following notation: 5kw is the diagonal matrix operator whose 
kth. diagonal entry equals multiplication by 5{y — C,kw), and 

kw 

It appears above because Dx = 6. 

^That it seems a miracle to us shows that we do not really understand why the equations should 
close. 
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Set 

p={I-Kx)-\ R = pK = {I ~ KXK)-^ - I. 
It follows from the last commutator upon left- and right-multiplication by p that 

[D, p] = -p(p(g)'iljxp + R6 p. 

From the commutators of D with K and p we compute 

[D, R] = [D,pK]=p [D, K] + [D, p]K = -py^^tp (/ + XR) + R6 R. 

Notice that 

I + XR= ii-xKy\ 

If we recall ()4.H] we see that we have shown 

[D, R] = -Q®P + R6R. (4.12) 

To obtain our second commutator we observe that if we apply dt + dg to the integrand 
in the formula for Hij{x, y) we get zero for the resulting integral. If we apply it to (s— 
we also get zero. Therefore we get zero if we apply it to the numerator, and this operation 
brings down the factor 

t"^ — Tit + X — S'^ + Tj s — y. 

The same factor results if we apply to Hij{x, y) the operator 

dl + dl-iT,d., + T,dy) + ix-y). 

We deduce 

[D^ -tD + M, H] = 0. 
One verifies that also [D^ — tD + M, E] = 0. Hence 

[D^ -tD + M, K] = 0. (4.13) 

This is the second commutator.^ From it we obtain 

[D^ -tD + M, Kx] = K [D^ -tD + M, X] = K {D6D + D^6 + 5D'^ - r5), 

and this gives 

[D'^ -tD + M, p] = R{D6D + D^6 + 6D^) p - Rt 6 p, (4.14) 

^ From 14.11|l we obtain also [D^, K] = —tp" ® if) -\- ip' ® ip' — ® i'" ■ Combining this with (|4.11|l 
itself and H4.13|l for m = 1 with r = ti we obtain [M, K] — p" ® ij; — Lp' ® ij;' + (p ® ij;" — ^ p ^ . This 
is equivalent to the formula stated in footnote [7| 
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which in turn gives 

[D^-tD+M, R] = [D^-tD + M, pK] = R{D5D + 0^5 + 50"^) R~Rt 6 R. (4.15) 



Of our nine equations the first seven are universal — they do not depend on the 
particulars of the kernel K or vector functions or tp. (The same was observed in |21j.) 
What are not universal are equations ()4.9|) and ()4.10p . For the equations to close we 
shall also have to show that the combinations of the entries of Vxx, r^y and Vyy which 
actually appear in the equations are all expressible in terms of the unknown functions. 
The reader can check that these are the diagonal entries of r^x + fxy and Vxy + Vyy (which 
also give the diagonal entries of Vxx — fyy) and the off-diagonal entries of Txx, ^xy and 

What we do at the beginning of our derivation will be a repetition of what was done 
in [21]. First, we have 

dkn,P = p{KdkwX)p={-lTR5kn,p. (4.16) 
From this we obtain dkwR = (—1)"' R^kw R, and so 

dkw ^iu,jv iS^kw Rij^ i,^iui ^jv') ~l~ ^xiji^^iui ^jv) ^iu,kw ~l~ Fiyij i,^iu^ ^j^^ ^jv,kw 
( 1) ^iu,kw ^kw,jv ~l~ ^x,iu,jv ^iu,kw ~l~ ^y,iu,jv ^jv,kw 

Multipliying by d^kw and summing over kw give ()4.2|1 . Equations ()4.Hj) and ()4.4j) are 
derived analogously. 

Next we derive ()4.5j) and ()4.7|1 . Using ()4.16|) applied to ip we obtain 

9kw Qiu Qii.^iu) ^iu,kw ~l~ ( 1) {^R^kw Q) ii^^iu) Qiu ^iu.kw ~l~ ( 1) ^iu,kw Qkw 

Multiplying by d^kw and summing over kw give ()4.5p . If we multiply ()4.16|) on the left 
by D we obtain dkwPx = {—^)^RxSkwP and applying the result to ip we obtain ()4.7p 
similarly. 

For ()4.9|) we begin as above, now applying to ()4.16|) on the left and ip on the right 
to obtain 

dkw Qiu Qi (.^iu) ^iu,kw ~l~ ( 1) '^xx,iu,kw Qkw (4'17) 

Now, though, we have to compute the first term on the right. To do it we apply ()4.14|) 
to ip and use the differential equation satisfied by (p to obtain 

Q"'{x) -tQ'{x) +xQ{x) = -Ry5Q' + Ryy5Q + R5Q" - Rt5Q. (4.18) 

This gives 

Qiiiiu) = Ti g.„ - Qiu + {-Ty sq' + ryysq + rsq" -rrs q)iu. 
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If we substitute this into ()4.17p . multiply by d^kw and sum over kw we obtain equation 
(gH). 

This completes the derivation of the equations for the differentials of g, q' and q" . 
We could say that the derivation of the equations for the differentials of p, p' and p" 
is analogous, which is true. But here is a better way. Observe that the P for the 
operator K is the transpose of the Q for the transpose of and similarly with P and 
Q interchanged. It follows from this that for any equation involving Q there is another 
one for P obtained by replacing K by its transpose (and so interchanging and dy) and 
taking transposes. The upshot is that equations ()4.6|) . ()4.8|) and ()4.10|) are consequences 
of ()4.5|) . ()4.7|) and ()4.9p . The reason for the difference in signs in the appearance of ^ 
on the right sides of ()4.9|1 and ()4.10|1 is the difference in signs in the last terms in the 
differential equations for (p and '0. 

Finally we have to show that the diagonal entries of Txx + ^xy and r^y + Vyy and 
the off-diagonal entries of r^x-, ^xy and Tyy are all known, in the sense that they are 
expressible in terms of the unknown functions. This is really the heart of the matter. 

We use = between expressions involving i?, Q and P and their derivatives to indicate 
that the difference involves at most two derivatives of Q or P and at most one derivative 
of R. The reason is that if we take the appropriate entries evaluated at the appropri- 
ate points we obtain a known quantity, i.e., one expressible in terms of the unknown 
functions. 

If we multiply fl4.12|l on the left or right by D we obtain 

Rxx ~l" Rxy ~ — Q' ® P ~\~ Rx ^ R^ Rxy ~l" ^yy ~ — Q ® P' ~\~ R S Ry. 

In particular 

Rxx ~t~ Rxy — 0; Rxy ~t~ Rxy — 0? 

SO in fact all entries of r^^ + r^.^ and r^y + Vyy are known. 
From ()4.12|) we obtain consecutively 

[D\ R] = -Q' 0P + Q0P' + DR6R + R6RD, (4.19) 
[D\ R] = -Q" 0P + Q' ®P' -Q0P" + D'^R6R + DR6RD + R6RD\ (4.20) 
If we subtract ()4.15j) from ()4.20j) we find 

[tD-M, R] = -Q" ®P + Q' ®P' -Q®P" 
-\-Rxx ^ R — Rx SRy -\- R6 Ryy + Ry 6 Rx — Ryy^ R — RSRxx -\- Rt 5 R. 

We use 

Rxx - Ryy = Q ® P' - Q' ® P + Rx5 R - R5 Ry 

and 

Ry = -R^-Q®P + RSR 

to see that this equals 
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-Q" P + Q' P' - Q P" + {Q P' - Q' P - R5 Ry) 5 R 

-R6{Q0P' -Q'0P + R^6R-R6Ry) + R^6Q0P (4.21) 
-{Q0P -R6R)6R^ + Rt6R. 

We first apply D acting on the left to this, and deduce that 

D [tD, R] = -Q'" 0P + R,,6Q0P. 

If we use fl4.18|l and the fact that Ryy = Rxx we see that this is = 0. This means that 

Since Rxx,i,j + Rxy,i,j = we deduce that Rxx,i,j and Rxy,i,j are individually = when 
i ^ j- Therefore rxx,iu,jv and r^y^i^jv are known then. 

But we still have to show that Viu^iv is known when u v, and for this we apply 
to (j4.2ip rather than D. We get this time 

D'^[tD~M, R] = -Q"" P + Q'" {P' - P 6 R) - Rxx6 Ry6 R 

- Rxx5 {Q P' - Q' P + R^S R - RS Ry) (4.22) 
+Rxxx 5Q0P + RxxSRSRx + RxxtSR. 

We first compare the diagonal entries of [t D, R] on the left with those of Rxx t 6 R 
on the right. The diagonal entries of the former are the same as those of 

T i^Rxxx ~l~ Rxxy ) = T Rxx ^ R- 

(Notice that applying to ()4.12j) on the left gives Rxxx + Rxxy = Rxx SR.) The 
difference between this and Rxx t6 Ris [r, Rxx] S R. Only the off-diagonal entries of Rxx 
occur here, so this is = 0. 

If we remove these terms from ()4.22|) the left side becomes —D^ [M, R] and the 
resulting right side we write, after using the fact Rx + Ry = —Q P + R6 R twice, as 

-Q"" 0P + Q'" 0{P' - P5R) 

- Rxx5 {Q P' - Q' ® P - Q ® P ^ R + Q ® P) + Rxxx5 Q P. (4.23) 

Now we use (j4.18|) and the facts 

Ryy ~ Rxx J Rxy — Rxxi Rxxx Rxyy — Rxx^ R 
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(the last comes from applying D to ()4.19|) on the left) to obtain 

Q'" = Rxx 5 Q, Q"" = Rxx S Q' + {Rxxx — Rxx 6 R) 6 Q. 

Substituting these expressions for Q'" and Q"" into ()4.23|) shows that it is = 0. 
This was the miracle. 
We have shown that 

[M, R] = 0, 

in other words (x — y) Rxx{x, y) = 0. If we set x = C,iu, y = C,iv we deduce that 
Rxx{^iu,^tv) = rxx,iu,iv is known when u y^v. 

V. Higher-order singularities 

We begin with the fictitious Brownian motion model, in which the end-points of the 
paths are complex numbers. The model consists of 2Rn nonintersecting Brownian paths 
starting at zero, with n of them ending at each of the points ±n^^'^br {r = 1, . . . , R). 
The product in the integrand in 1)2.111) becomes 

and we use the same contours as before. 

We shall first make the substitutions 1/2 + n^^rfc with 6 to be determined. The 
first exponential in ()2.11|) becomes 

exp { - (1 + 4:n-^Tk + 0{n-^^)) + (4 + 0{n-^)) xt + ©(x^)}, (5.2) 

and the second exponential becomes 

exp {(1 + 4n-^n + 0{n-^^)) s^-{A + 0{n-^)) ys + 0{y^)]. 

If we set Oj. = 1/6^ the product ()5.1|) is the exponential of 



1 "^^^ 

+ 



+0(n-«-^(|t|^^+^ + |s|^«+^)). 
If i? > 1 we choose the such that 

^0^ = 1, = ■ ■ ■ = = 0. 
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The ttr are the roots of the equation 

- a^-' + ia^-^ - ■ ■ ■ + = o, (5.3) 

from which it follows that 

In general the will be complex, and so the same will be true of the end-points of our 
Brownian paths. 

In the integrals defining we make the substitutions 



where 



R + r 

and in the kernel we make the scaling 

X rT^^"^ X, y rT^^"^ y. 

This gives us the kernel Hn,ij{x,y). As before Eij is unchanged, and the limiting form 
of Hn,ij{x,y) is now 



g(-l)^s2«+2/{iJ+l)!+4T,- s'^-Ays+{-l)^+^t^^+^/{R+iy.-ATit^+4xt "^"^ 



Jc J-ioo S — t 

This is formal and it is not at first clear what the C contour should be, although one 
might guess that it consists of four rays, one in each quadrant, on which (_i)-f?+it2R+2 
is negative and real. We shall see that this is so, and that the rays are the most vertical 

I _H_ ■ 

ones, those between and icxoe 2^+2 . xhe orientation of the rays is as in the case 
R= 1. (The s integration should cause no new problems.) 

After the variable changes the product of the two functions of t in the integrand in 
()2.11|) is of the form 



The main part of this is the quotient. 

Upon the substitution t n^^^^^^'^t the quotient becomes 



(5.5) 
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Suppose we want to do a steepest descent analysis of the integral of this over a nearly 
vertical ray from in the right half-plane. (This nearly vertical ray would be the part 
of the t contour in ()2.1Hl in the first quadrant.) No pole ±6^ is purely imaginary, as 
is clear from the equation the satisfy. So there are R poles in the right half plane 
and R in the left. There are 2R + 2 steepest descent curves emanating from the origin, 
half starting out in the right half-plane. These remain there since, as one can show, 
the integrand is positive and increasing on the imaginary axis. We claim that there is 
at least one pole between any two of these curves. The reason is that otherwise the 
integrals over these two curves would be equal, and so have equal asymptotics. That 
means, after computing the asymptotics, that the integrals 

ifc7r/(R+l) 
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would be the same for two different integers k G (— (-R + l)/2, {R + l)/2). But the 
integrals are all different. 

Therefore there is a pole between any two of the curves. Let F be the curve which 
starts out most steeply, in the direction arg t = ^^^ tt. It follows from what we have 
just shown that there is no pole between T and the positive imaginary axis. This is what 
we wanted to show. 

The curve we take for the t integral in ()2.11|) is r„ = n^^^^y^ T. The original contour 
for the t-integral in the representation of Hn,ij can be deformed to this one. (We are 
speaking now, of course, of one quarter of the full contour.) 

We can now take care of the annoying part of the argument establishing the claimed 
asymptotics. The curve F is asymptotic to the positive real axis at +oo. Therefore for 
A sufficiently large is 0(e-"l*l'/2) when t G T, \t\ > A. Hence dEl is 0(e-"'l*l'/^) 
when t G r„, \t\ > n^^'^^^^A. It follows that its norm over this portion of r„ is 
exponentially small. When t G F, |t| > e then ()5.5p is 0(e~"'') for some rj > 0, and 
it follows that (Q is 0(e-"''/2) ^^g^ ^ ^ y^, \t\ > d^'^'>/^e and also \t\ < d^'^^/'^A. 
Therefore the norm of (j5.4p over this portion of r„ is also exponentially small. So we 
need consider only the portion of r„ on which |t| < n^^~^^^'^e, and for this we get the 
limit 

Jo 

with appropriate uniformity, in the usual way. 

Just as with the Pearcey kernel we can search for a system of PDEs associated with 
det {I — Kx). Again we obtain two commutators, which when combined show that K 
is an integrable kernel when m = 1. In this case we define (p and ip by 



Kni Jc 
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Xiri J-ioo J 
They satisfy the differential equations 

where 

The first commutator is [D, K] = —4^9 (g) ip, which also gives 

ra-l 

[D", ir] = -4 ^ (-f ^(""'=-1) ® ip^'^l (5.6) 

fe=0 

The second commutator is 

[cr D^^+^ -2tD + AM, K]=0. 

In case m = 1 (or for a general m and a diagonal entry of K), combining this with the 
commutator [D"^^^^, K] obtained from ()5.(jj) and the differential equations for ip and ip 
we get an expression for [M, K] in terms of derivatives of (p and ijj up to order 2R. This 
gives the analogue of the expression for K{x, y) in footnote H 

For a system of PDEs in this case we would have many more unknowns, and the 
industrious reader could write them down. However there will remain the problem of 
showing that certain quantities involving 2i?th derivatives of the resolvent kernel R (too 
many i?s!) evaluated at endpoints of the Xk are expressible in terms of the unknowns. 
For the case i? = 1 a miracle took place. Even to determine what miracle has to take 
place for general R would be a nontrivial computational task. 
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